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Abstract 

A general model of catalytic branching process (CBP) with any finite number of catal¬ 
ysis centers in a discrete space is studied. More exactly, it is assumed that particles move 
in this space according to a specified Markov chain and they may produce offspring only 
in the presence of catalysts located at fixed points. The asymptotic (in time) behavior of 
the total number of particles as well as the local particles numbers is investigated. The 
problems of finding the global extinction probability and local extinction probability are 
solved. Necessary and sufficient conditions are established for phase of pure global survival 
and strong local survival. Under wide conditions the limit theorems for the normalized 
total and local particles numbers in supercritical CBP are proved in the sense of almost 
surely convergence as well as with respect to convergence in distribution. Generalizations 
of a number of previous results are obtained as well. In the proofs the main role is played 
by recent results by the author devoted to classification of CBP and the moment analysis 
of the total and local particles numbers in CBP. 
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Theory of branching processes is a classical part of probability theory (see, e.g., H)- Its 
origin goes back to the middle of XIX centnry when F.Galton and G.W.Watson proposed a 
model elncidating the extinction of family names. In the framework of this model the pri¬ 
mary branching process called Galton-Watson process arose. It was fonnd that the valne of 
extinction probability as well as the mean popnlation size depend essentially on whether the 
mean offspring nnmber of a popnlation representative is greater than, eqnal to or less than 1. 
Galton-Watson branching process having that mean offspring nnmber is called snpercritical, 
critical and snbcritical, respectively. Only for snpercritical Galton-Watson process the extinc¬ 
tion probability is less than 1 whereas the mean popnlation size tends to inhnity as time grows. 
Therefore, after solving the problem of Ending the extinction probability the qnestion arises 
how fast the popnlation grows in the case of its snrvival. The answer to that qnestion is given 
by limit theorems for the popnlation size. The resnlts related to the Galton-Watson process 
are contained, e.g., in la, sections 1-6 and 9. 

We are interested in a more involved model than the Galton-Watson process, namely cat¬ 
alytic branching process (GBP) with an arbitrary finite nnmber of catalysis centers. Its specifics 
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consists in possibility for the population representatives (which further will be called particles) 
not only to produce offsprings but also to move in the space. Moreover, we assume that the par¬ 
ticles produce offsprings exclusively in the presence of catalysts which are located in a hnite set 
of the space points. For CBP it is natural to raise the question not only on the global particles 
extinction but also on the local extinction as well as to investigate the asymptotic behavior of 
total and local particles numbers. The present paper is devoted to these problems. We consider 
a discrete space where the particles move, although branching processes in a continuous space 
are of doubtless interest as well (see, e.g., 0). 

The study of general model of CBP with any finite number of catalysis centers in a discrete 
space is initiated in |3]. Particular cases of this model were analyzed in many publications, 
e.g., in n, 13,0, HI and others. In dependence on the value of the Perron root of a certain 
matrix a classihcation of CBP as supercritical, critical and subcritical processes was proposed 
in |1]. The naturalness of the classification was confirmed by implemented by the author 
asymptotic analysis of the moments of the total and local particles numbers existing in the 
process at time t, as t —?• cxo. In the present work the problem of finding the probabilities of 
global and local population extinction is solved and new limit theorems in strong (Theorem 
|3]) and weak (Theorem 0]) forms are established for the total and local particles numbers in a 
supercritical CBP. The results obtained generalize a number of previous ones proved, e.g., in 
papers [7] (Lemma 5.1) and [U]. So, limit distributions for particles numbers in the model 
of supercritical branching random walk on integer lattice Z'^, d G N, with a finite number 
of particles generation centers are found in |T5]. There not only existence of all moments of 
an offsprings number for any particle was assumed but also a certain rate of growth of such 
moments in dependence on its order. In our work, besides considering a more general model 
we also impose extremely weak restrictions on the moments of the offsprings size for each 
particle and, moreover, study the limit behavior of the particles numbers in the sense of almost 
surely convergence. These results are not only of self-contained interest but will be also applied 
in further investigations of the character of the particles population propagation in CBP. It 
is worthwhile to note that the mentioned achievements in the domain of limit theorems for 
supercritical CBP are due to the foundation laid in a recent work [1]. Firstly, there were 
introduced auxiliary Bellman-Harris branching processes which mainly allows us to reduce the 
study of CBP to classic results of modern theory of branching processes. Secondly, the moment 
analysis of the total and local particles numbers was implemented in |1] and the proofs of results 
of the present paper substantially bear on it. 

Let us give a formal description of CBP. At the initial time t = 0 there is a single particle, 
its movement on a specihed hnite or countable set S is governed by continuous-time Markov 
chain rj = {r](t),t > 0} with inhnitesimal matrix Q = {q{x,y))x,y£S- Hitting a hnite catalysts 
set W = {wi ,... ,wn} C S, e.g., at point Wk, the particle spends there random time distributed 
according to the exponential law with parameter 13k > 0. Then it either produces ohsprings 
or leaves the element Wk with probabilities and 1 — (0 < < 1), respectively. If 

the particle produces ohsprings (at Wk), it dies instantly turning into a random number ^k of 
ohsprings located at the same point Wk- If the particle leaves Wk, it jumps to point y ^ Wk 
with probability —q{wk,y)q{wk,Wk)~^ and continues the movement controlled by the Markov 
chain y. All the newborn particles behave themselves as independent copies of the parent 
particle. We assume that the Markov chain r] is irreducible and the matrix Q is conservative, 
i.e. = 0 as well as q{x,y) > 0 for x ^ y and q{x,x) G (— cxd,0) for any a; G S'. 

Denote by fk{s) := Es^^, s G [0,1], the probability generating function of random variable ^k, 
k = 1,..., A^. We will employ the standard condition of the hnite derivative /(.(I) existence. 
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i.e. of finiteness for each k = 1,..., N. Let fi{t) be the total particles number existing in 
CBP at time f > 0 and the local particles numbers /i(f; y) are the quantities of particles located 
at separate points y E S a.t time t. 

To formulate the main results of the paper let us introduce additional notation. We tem¬ 
porarily forget that there are catalysts at some elements of set S and consider only the movement 
of a particle on set S in accordance with Markov chain r] with generator Q and starting point 
X. Let HTx,y be the time elapsed from the exit moment of this Markov chain out of starting 
state X till the first hitting point y by the particle whenever the particle trajectory does not 
pass the set H <Z S. Otherwise, we put HTx,y = oo. Extended random variable Hfx,y is called 
hitting time of state y under taboo on set H after exit out of starting state x (see, e.g., [8], 
Part II, Section 11). Denote by HFx,y{t), t > 0, the improper cumulative distribution function 
of this extended random variable. Mainly we will be interested in the situation when H = Wk 
where Wk := W \ {wk}, k = 1,..., N. Further 

poo 

F*{X):= / A > 0, 

Jo- 


denotes the Laplace transform of cumulative distribution function F{t), t > 0, with support 
located on non-negative semi-axis. For j,k = 1,... ,N and f > 0 set Gj{t) := 1 — 



Gj^kijJ) ■— /Jj j WkFwj,wki^ u)e du 

and 

WkFx,wk{t) :=-q{x,x) [ du 

Jo 

At last, let 


WkFx,Wk{^) — WkFx,Wk{^) ■— fil^ WkFx,Wk{f)- 

t^OO 


In [4] there was introduced a matrix function D{X) which is irreducible matrix of size N x N 
for each A > 0. Namely, D(A) = {dij{X))F^i where 


dij{X) — 6ijaifl{l)G*{X) + (1 - ai)G*{X)wjF.^. .^^{X) 


and 6ij is the Kronecker delta. According to definition 1 in [1] CBP is called supercritical if 
the Perron root (i.e. positive eigenvalue being the spectral radius) p{D{0)) of the matrix D(0) 
is greater than 1. Then in view of monotonicity of all elements of matrix function D{-) there 
exists the solution u > 0 of equation p{D{X)) = 1. As Theorem 1 in [1] shows, just this positive 
number z/ specifies the rate of exponential growth of the mean total and local particles numbers 
(in the literature devoted to population dynamics and branching processes one traditionally 
speaks of Malthusian parameter). More exactly, Exp(t) ~ A{x)e’'^ and Exp(t;y) ~ a{x,y)e’^^ as 
t ^ oo where the index x E S stands for the starting point of CBP. The explicit formulae for 
functions A(-) and a(-, •) are given in [1]. 

Theorem [1] is devoted to solution of the problem of finding the global extinction probability 
q{x) ;= Pa,(hmt_,.oo /i(t) = 0) = hmi_j.oo Pxihit) = d), x E S, for the model of CBP. 


Theorem 1 For the global extinction probability q{x) when x E S\W there exists a represen¬ 
tation 

N 

= '^w^Fx,,u,^{oo)q{wk) ( 1 ) 

k=l 
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where the values q{wj), j = 1,..., satisfy with the following system of equations 

N 

q{wj) = ajfj{q{wj)) + (1 - aj) ^ w^F^j,w^{oo)q{wk). ( 2 ) 

k=l 

In addition, vector {q{wi),..., q{wN)) is component-wise the least root of equations system ([2]) 
in the cube [0,1]^. Moreover, if the Markov chain rj is recurrent then q{x) = 1, x E S, or 
q{x) < 1, X E S, whenever p{D{0)) < 1 or p{D{0)) > 1, respectively. If the Markov chain p is 
transient then q{x) < 1 for all x E S. 

One may oppose the global extinction probability q{x) to the local extinction probability 
Q{x,y) = Pa; (limsnp^^oo/i(t;?/) = 0), for x,y E S. The following theorem shows that in fact 
the fnnction Q{x,y) does not depend on variable y. 

Theorem 2 Equality Q{x,y) = Q{x) holds for any y E S where the function Q(x) when 
X E S \ W is of the form 


N 


N 


X,W}^ 


.00 


Q(^x') ^ ^ (oo)Q(^fc) T 1 ^ ^ 

k=l k=l 

and the values Q{wj), j = 1,..., N, are the least solution to the equations system 

N 

Qi'^j) ~ ^jfji.Qi.'^j)) T (1 — dj) ^ ^^ WkFwj,wk{'^)Q{'^k) 

k=l 

/ N '' 

+ (1 — otj) 11 — 'y ^ WkFwj,wk{’^) 


k=l 


( 3 ) 


( 4 ) 


in the cube [0,1]^. Moreover, if p{D{0)) < 1 then Q{x) = 1 for all x E S. Whenever 
p{D{0)) > 1 one has Q{x) < 1 for each x E S. 


Applying Theorems 4 and 5 in [12], Chapter 5, Section 1, to the anxiliary Bellman-Harris 
processes constructed in [1] we come to the statements of Theorems [T] and |2l Clearly 0 < 
q{x) < Q{x) < 1, X E S. In view of the explicit form of relations (11])-(Hj) we conclude that if 
the strict inequality '^^=1 WkFx,wk{'^) < 1 is satished at least for some x E S, i.e. the Markov 
chain p is transient, then q{x) < Q{x) for all x E S. Otherwise, i.e. if the Markov chain p is 
recurrent, the relations for g(-) and Q{-) coincide and whence q{x) = Q{x) for all x. In the 
terms of paper |2] the aforesaid means that, for transient Markov chain p and p{D{0)) < 1, we 
deal with the pure global survival phase of CBP and in the case of recurrent Markov chain p 
and p(Zl(0)) > 1 the strong local survival of CBP is observed. 

Now we pass to considering the problem of the population growth rate in the case of global 
and local survival, i.e. whenever p{D{0)) > 1. Let u = (ui,... ,un) be the right eigenvector 
of matrix D{i>) corresponding to the Perron root p{D{v)) equal to 1 where > 0, fc = 
1,..., A^, and 'Wfc = 1- If should be noted that, by virtue of the Perron-Frobenius theorem 
(see, e.g., [12], Chapter IV, Section 5), such eigenvector can always be found since the matrix 
D{v) is irreducible according to Lemma 1 in [1]. Recall the dehnition of matrix function 
D{x]\), X ^ W, A > 0, introduced in [1] while proving Theorem 1, case 2. For x ^ W, set 
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wn +1 = X, W{x) ■.= W U {x} and Wi{x) := W{x) \ {wi}, i = 1,..., N + 1. Then the matrix 
D{x; A) = {dij{x; has elements 


A) := 5,, aj:{l)G*{X) + (1 - a,)G:{X)w,i.)Fl(X), X > 0. 


Here oat+i = 0, = 0 and GAr+i(t) := 1 — t > 0. By Lemma 3 in |1] the matrix 

D{x-, A) is irredncible and p{D{x] v)) = 1. Define u(x) = (mi(x), ... ,MAr+i(x)) to be the right 
eigenvector of the matrix D{x] z/), corresponding to the Perron root p{D{x] v)) eqnal to 1, snch 
that Uk{x) >0, A; = 1,. .., iV + 1, and Ylk=i '^k{x) = 1. Take notice that bearing on the proof 
of Lemma 3 in [1] it is not difficult to verify equalities Ui{x) = Ui for all i = 1,... ,iV and 
X E S \ W. Set c(x) := whenever x = Wk for some k = 1,.. .,N and c(x) := uX/j^_^{x) 

whenever x G S'\iy. Let also symbols 0 and 1 denote vectors with zeros and units, respectively, 
as all components, the dimension of vectors being contextually clear. 

In the given paper, we handle three forms of random variables convergence, viz., almost 
surely, in probability and in distribution which are denoted by A- and A, respectively. 

Theorem [3] describes strong convergence of vectors of the normalized total and local particles 
numbers in CBP as time t grows to infinity. The results on asymptotic behavior of the normal¬ 
izing means obtained in Theorem 1 in |3] are implicitly employed in the proofs of Theorems [3] 
andm 


Theorem 3 Let a supercritical CBP start at point x E S. Assume that the elements of matrix 
Q are uniformly bounded, i.e. forallzi,Z 2 E S and some constantG > 0 one has \q{zi, Z 2 )\ < G. 
Let also < oo for each k = 1,..., N. Then there exists non-degenerate random variable ( 
such that for any n G M and yi,... ,yn E S the following relation is valid 


p{t) p(t,yi) p{t,yn) A a.s. , 

Ea;/r(f) ’ i/i) ’ ■ ■ ■ ’ Ea,/i(f; |/„) / 


t ^ oo. 


(5) 


If the Markov chain p is recurrent, the assertion of Theorem [3] ensues from the construction 
of auxiliary Bellman-Harris branching process in |1] and applying to it Theorem 4.1 in m with 
due regard of the fact that under the conditions of Theorem [3] the distributions of life-lengths 
of particles of different types in the auxiliary Bellman-Harris process are absolutely continuous 
and have bounded densities. The latter observation is true on account of the argument in the 
final part of paper [5]. 

If the Markov chain p is transient then one has to apply the proof scheme of Theorem 4.1 in 
m to the auxiliary Bellman-Harris process with final type of particles which was introduced 
in [1] while establishing the case 1 (step 6) of Theorem 1. 

Now we relax the restrictions on the studied process. For 2 ; G M, set log"*" 2 ; := log(max{^, 1}). 


Theorem 4 Let a supercritical CBP start at point x E S. 
yi, ■ ■ ■ ,yn E S the following alternative is true. 

1. If E^k log^ fk = oo for some G {1,..., N} then 

/ /x(t) /x(t;^i) ^ 

yn) J 


Then for each n E N and any 



t ^ 00 . 


( 6 ) 


2. If E^k log’*' fk<oo for all k = 1,..., N then 

f T{t) Likyi) L{t-,yn) A 

V E„/r(t) ’ E,,p{t] 2 / 1 ) ’ ■ " ’ E,,p{t] yn) J 


4 c(x)Cl, 


t ^ 00 . 


(7) 
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Here ( is a non-degenerate random variable with the following properties. 

(i) E^C = cix)-^. 


(ii) P^(C = 0) = (limsupi^^/x(t;?/) = 0) for any y E S. 

(iii) The Laplace transform ip{X; x) := \ > t), x E S, of random variable ( forx E S'\PP 

is of the form 


^ POO 

(p{X]x) = ^ / (p{Xe~'''^]Wk) dwuFx,wk{'^) + 1 

fc=i 


N 

^ ^ VKfc-^a;,«)fe(oo) 

k=l 


where functions ip{--,Wj), j = 1,... ,N, satisfy the system of integral equations 

POO 

^{X]Wj) = aj / fj{^p{Xe~'^'';Wj))dGj{u) 

Jo 

^ POO 

+ (1 - aj) ^ / (p(Ae“''“; Wk) dGj^k{u) 


k=l 


N 


+ (1 Oj) j 1 ^^WkFwj,vn.{oo) 


k=l 


(iv) The conditional distribution of C, under condition of CBP start at point x E S is absolutely 
continuous on the positive semi-axis and has continuous density function. 

Note that in relation ((71) of Theorem 0] exactly the convergence of vectors is essential whereas 
in formulae ([5]) and ((6]) the vectors convergence is tantamount to convergence of their compo¬ 
nents. 

If the Markov chain r] is recurrent then Theorem 0] is proven with the help of Theorem 1.1 in 
[To] , applied to the constructed in |1] auxiliary Bellman-Harris process. However, if the Markov 
chain p is transient, the statement of Theorem 0] is established by applying the proof scheme of 
Theorem 1.1 in ra to the Bellman-Harris process with hnal type of particles introduced while 
proving the case 1 (step 6) of Theorem 1 in [1]. 

Thus, under weak conditions new asymptotic properties of CBP are investigated. 
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